Coulomb-assisted braiding of Majorana fermions in a Josephson junction array 
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We show how to exchange (braid) Majorana fermions in a network of superconducting nanowires 
by control over Coulomb interactions rather than tunneling. Even though Majorana fermions are 
charge-neutral quasiparticles (equal to their own antiparticle), they have an effective long-range 
interaction through the even-odd electron number dependence of the superconducting ground state. 
The flux through a split Josephson junction controls this interaction via the ratio of Josephson and 
charging energies, with exponential sensitivity. By switching the interaction on and off in neighboring 
segments of a Josephson junction array, the non-Abelian braiding statistics can be realized without 
the need to control tunnel couplings by gate electrodes. This is a solution to the problem how to 
operate on topological qubits when gate voltages are screened by the superconductor. 

PACS numbers: 03.67.Lx, 73.23. Hk, 74.50.+r, 74.81.Fa 



I. INTRODUCTION 

Non-Abelian anyons have a topological charge that 
provides a nonlocal encoding of quantum information pQ . 
In superconducting implementations [2] |3] the topologi- 
cal charge equals the electrical charge modulo 2e, shared 
nonlocally by a pair of midgap states called Majorana 
fermions [3]. This mundane identification of topologi- 
cal and electrical charge by no means diminishes the rel- 
evance for quantum computation. To the contrary, it 
provides a powerful way to manipulate the topological 
charge through the well-established sub-e charge sensi- 
tivity of superconducting electronics [6] . 

Following this line of thought, three of us recently 
proposed a hybrid device called a top-transmon, which 
combines the adjustable charge sensitivity of a super- 
conducting charge qubit (the transmon [7]) to read out 
and rotate a topological {top) qubit [8.. A universal 
quantum computer with highly favorable error thresh- 
old can be constructed [5] if these operations are supple- 
mented by the pairwise exchange (braiding) of Majorana 
fermions, which is a non-Abelian operation on the degen- 
erate ground state (TUJ [IT] . 

Here we show how Majorana fermions can be braided 
by means of charge-sensitive superconducting electron- 
ics. (Braiding was not implemented in Ref. [8 nor in 
other studies of hybrid topological/nontopological super- 
conducting qubits P^rfTB] .) We exploit the fact that the 
charge-sensitivity can be switched on and off with expo- 
nential accuracy by varying the magnetic flux through a 
split Josephson junction |7J. This provides a macroscopic 
handle on the Coulomb interaction of pairs of Majorana 
fermions, which makes it possible to transport and ex- 
change them in a Josephson junction array. 

We compare and contrast our approach with that of 
Sau, Clarke, and Tewari, who showed (building on the 
work of Alicea et al. [T7]) how non-Abelian braiding 
statistics could be generated by switching on and off the 
tunnel coupling of adjacent pairs of Majorana fermions 
[18] , The tunnel coupling is controlled by a gate volt- 



age, while we rely on Coulomb interaction controlled by a 
magnetic flux. This becomes an essential difference when 
electric fields are screened too strongly by the supercon- 
ductor to be effective. (For an alternative non-electrical 
approach to braiding, see Ref. [TTj].) 

The basic procedure can be explained quite simply, see 
Sec. |III| after the mechanism of the Coulomb coupling is 
presented in Sec. |H] We make use of two more involved 
pieces of theoretical analysis, one is the derivation of the 
low-energy Hamiltonian of the Coulomb coupled Majo- 
rana fermions (using results from Refs. [20][2T]), and the 
other is the calculation of the non-Abelian Berry phase 
[22] of the exchange operation. To streamline the paper 
the details of these two calculations are given in Appen- 
dices. 



II. MAJORANA-COULOMB HAMILTONIAN 

A. Single island 

The basic building block of the Josephson junction ar- 
ray is the Cooper pair box [23], see Fig. [Tl consisting 
of a superconducting island (capacitance U) connected 
to a bulk (grounded) superconductor by a split Joseph- 
son junction enclosing a magnetic flux $. The Joseph- 
son energy Ej is a periodic function of $ with period 
$o = h/2e. If the two arms of the split junction are 
balanced, each with the same coupling energy Eq, the 
Josephson energy 

Ej = 2E cos(7T$/$ ) (1) 

varies between and 2Eq > as a function of |<E>| < $0/2. 

When the island contains no Majorana fermions, its 
Hamiltonian has the usual form [2JJ 

H= ^(Q + qind) 2 -Ej cos 0, (2) 

in terms of the canonically conjugate phase <fi and charge 
Q = —2eid/d(f> of the island. The offset (find accounts 
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FIG. 1: Cooper pair box, consisting of a superconducting 
island (brown) connected to a bulk superconductor by a split 
Josephson junction (black, with the gauge-variant phase dif- 
ferences indicated). The island contains Majorana fermions 
(yellow) at the end points of a nanowire (grey). These are 
coupled by the Coulomb charging energy, tunable via the flux 
<f> through the Josephson junction. 

for charges on nearby gate electrodes. We have chosen a 
gauge such that the phase of the pair potential is zero on 
the bulk superconductor. 

A segment of a semiconductor nanowire (typically 
InAs) on the superconducting island can have Majorana 
midgap states bound to the end points [21 [3] . For N seg- 
ments there can be 2N Majorana fermions on the island. 
They have identical creation and annihilation operators 
In = it satisfying 

7" 7m + Tm7n = 2£ n m< (3) 

The topological charge of the island equals the fermion 
parity 

2JV 

V = i N J] 7»- (4) 

n=l 

The eigenvalues of V are ±1, depending on whether there 
is an even or an odd number of electrons on the island. 

The Majorana operators do not enter explicitly in H, 
but affect the spectrum through a constraint on the 
eigenstates [20] . 

+ 2tt) = (-l)^- v ^H((f>). (5) 

This ensures that the eigenvalues of Q are even multiples 
of e for V = 1 and odd multiples for V = — 1. Since V 
contains the product of all the Majorana operators on 
the island, the constraint ^ effectively couples distant 
Majorana fermions — without requiring any overlap of 
wave functions. 

We operate the Cooper pair box in the regime that 
the Josephson energy Ej is large compared to the single- 
electron charging energy Ec — e 2 /2C. The phase 
<p (modulo 2tt) then has small zero-point fluctuations 
around the value m in — which minimizes the energy 
of the Josephson junction, with occasional 2tt quantum 
phase slips. 




FIG. 2: Two Cooper pair boxes, each containing a pair of 
Majorana fermions. Single electrons can tunnel between the 
superconducting islands via the overlapping Majorana's 712 
and 721 • This tunnel coupling has a slow (cosine) dependence 
on the enclosed fluxes, while the Coulomb coupling between 
the Majorana's on the same island varies rapidly (exponen- 
tially). 

In Appendix [A] we derive the effective low-energy 
Hamiltonian for Ej ^> Ec, 

H cS = -Ej + y/2E c Ej - UV, (6) 
U = miEcEj^n^e-V^j/Ea cos (^. nd / e ). (7) 

The energy minimum —2Eq at 4> m i n is increased by 
\J2EqEj due to zero-point fluctuations of the phase. 
This offset does not contain the Majorana operators, so 
it can be ignored. The term — UV due to quantum phase 
slips depends on the Majorana operators through the 
fermion parity. This term acquires a dynamics for mul- 
tiple coupled islands, because then the fermion parity of 
each individual island is no longer conserved. 

B. Multiple islands 

We generalize the description to multiple supercon- 
ducting islands, labeled k — 1,2,..., each connected to 
a bulk superconductor by a split Josephson junction en- 
closing a flux (See Fig. [2]) The Josephson junctions 
contribute an energy 

Hj = - ^ E J>k cos fa, E Jtk = 2E cos(7T$ fc /$ )- (8) 

k 

We assume that the charging energy is dominated by 
the self-capacitance C of each island, so that it has the 
additive form 

Hc=Y,^c^+1^) 2 - ( 9 ) 

k 

While both Eq and C may be different for different is- 
lands, we omit a possible fc-dependence for ease of no- 
tation. There may be additional fluxes enclosed by the 
regions between the islands, but we do not include them 
to simplify the expressions. None of these simplifications 
is essential for the operation of the device. 

The set of Majorana's on the fc-th island is indicated 
by 7fe n with n = 1,2, ...2N k . The fermion parities 
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FIG. 3: Three Cooper pair boxes connected at a tri-junction 
via three overlapping Majorana fermions (which effectively 
produce a single zero-mode 70 at the center). This is the 
minimal setup required for the braiding of a pair of Majo- 
rana's, controlled by the fluxes through the three Josephson 
junctions to a bulk superconductor. 



Yin 7 fe ™ °f neighboring islands k and k! are cou- 



pled with strength Em by the overlapping Majorana's 
7/c n and Jk'm- We denote the gauge-invariant phase dif- 
ference [21] by kk , = 4> k - <t>k> + (2tt/$ ) J k ^ k , A ■ dl. 
The corresponding tunnel Hamiltonian [3] 



H k k' = r fcfe ' cos(0 fc fc//2), r fcfc ' = iEMlknlk 1 ' 



(10) 



is 47r-periodic in the gauge-invariant phase difference, as 
an expression of the fact that single electrons (rather than 
Cooper pairs) tunnel through the midgap state. For ex- 
ample, in the two-island geometry of Fig. [2] one has 



H 12 = iE M li2l2i cos(6»i 2 /2), 

\ - 4> 2 - 7r($i + $ 2 )/$0- 



'12 



(11a) 
(lib) 

In Appendix [X] we derive the effective low-energy 
Hamiltonian in the regime Ej Eq, Em, 

H cff = const - ^2 u kVk + ^2 Tkk ' cos akk ' ' ( 12 ) 



a kk , = lim \Q kk >. 



0fc,0fc'->-O 



(13) 



The single sum couples Majorana's within an island, 
through an effective Coulomb energy U k . The double 
sum couples Majorana's in neighboring islands by tun- 
neling. Both the Coulomb and tunnel couplings depend 
on the fluxes through the Josephson junctions, but in an 
entirely different way: the tunnel coupling varies slowly 
ex cos(7T$/$ ) with the flux, while the Coulomb coupling 
varies rapidly oc exp[— Ay/ (E /E c ) cos(7r$/$o)]. 



will demonstrate, it can provide the non-Abelian braiding 
statistic that we are seeking. 

Similarly to earlier braiding proposals [171118] . the min- 
imal setup consists of three superconductors in a tri- 
junction. (See Fig. [3]) Each superconductor contains 
a pair of Majorana fermions 7fc,7j,, with a tunnel cou- 
pling betwee n 7^ , 7 2 , and 7 3 . The Majorana- Coulomb 
Hamiltonian |l2| takes the form 

H cS = iE M (7W2 cos "12 + 7 2 73 cos "23 + 7 3 7i cos a 3 i) 
3 

-$> fc *7fc7fc, (14) 
fc=i 

with gauge-invariant phase differences 



ai2 = -(tt/2$ )(*i + *2 
a 23 = (7r/2$o)($2 + $3), 
a 31 = (7r/2$ )($i + $ 3 ). 



2$ 3 ), (15a) 
(15b) 
(15c) 

As we vary \& k \ between and $ max < 3>o/2, the 
Coulomb coupling U k varies between two (possibly k- 
dependent) values ?7 m i n and t/ max . We require U max ^> 
U m in, which is readily achievable because of the exponen- 
tial flux sensitivity of the Coulomb coupling expressed 
by Eqs. ([I]) and We call the Coulomb couplings 

[/„,„ and {7 m i n on and off, respectively. We also take 
[/„,„ <C Em, meaning that the Coulomb coupling is 
weaker than the tunnel coupling. This is not an essential 
assumption, but it allows us to reduce the 6-Majorana 
problem to a 4-Majorana problem, as we will now show. 

Consider first the case that U k — for all k. Then 



the Hamiltonian ( 14 1 has four eigenvalues equal to zero: 



three of these represent the Majorana's j k far away from 
the junction, while the fourth Majorana, 

(16) 



7o 



1 W+72 +73) 



V3 



is situated at the tri-junction. The tri-junction con- 
tributes also two nonzero eigenvalues ±i-E gap , separated 
by the gap 



E gap = E M ycos- < i u - < •( >s- n 2 . ( 



COS z CX31. 



(17) 



For $ max well below $ and U mseK <C E M these two 
gapped modes can be ignored, and only the four Ma- 
jorana's 70,71,72,73 need to be retained. 

The Hamiltonian H mt that describes the Coulomb in- 
teraction of these four Majorana's for nonzero U k is given, 
to first order in U k /EM, by 

3 

H mt =J2 Ak hoik, A fe = -{2E M / Eg ap )f3 k U k , (18) 



fc=i 



/?i=cosa 23 , /3 2 =cosa 3 i, /3 3 = cosai 2 . 



(19) 



C. Tri-junction 



III. MAJORANA BRAIDING 



Since Vk and T kk > in the Majorana- Coulomb Hamil- 



tonian (12 1 do not commute, the evolution of the eigen- 
states upon variation of the fluxes is nontrivial. As we 



The Hamiltonian ( 18 1 describes four flux-tunable 



Coulomb-coupled Majorana fermions. Although the cou- 
pling studied by Sau, Clarke, and Tewari [18j has an 
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FIG. 4: Schematic of the three steps of the braiding op- 
eration. The four Majorana's of the tri-junction in Fig. [3] 
(the three outer Majorana's 71 , 72 , 73 and the effective cen- 
tral Majorana 70) are represented by circles and the Coulomb 
coupling is represented by lines (solid in the on state, dashed 
in the off state). White circles indicate Majorana's with a 
large Coulomb splitting, colored circles those with a vanish- 
ingly small Coulomb splitting. The small diagram above each 
arrow shows an intermediate stage, with one Majorana de- 
localized over three coupled sites. The three steps together 
exchange the Majorana's 1 and 2, which is a non-Abelian 
braiding operation. 



time 


$1 


$ 2 


$3 











-$max 




^max 





-$max 


T 


$max 










^max 


^*max 





2T 





^raax 










^*max 


-$max 


3T 








-*max 



TABLE I: Variation of the flux through the three Josephson 
junctions during the braiding operation, at time steps corre- 
sponding to the diagrams in Fig. [4] The flux $3 is varied in 
the opposite direction as $1, $2, to ensure that the coupling 
parameters oc /3k do not change sign during the operation. 



braiding operation then consists, in sequence, of the three 
steps 03i, 0i2, and 2 3. Note that each coupling 
appears twice in the on state during the entire operation, 
both times with the same sign Sfc. 

The step Okk' transfers the uncoupled Majorana at 
site k' to site k in a time T. The transfer is described in 
the Heisenberg representation by 7k (T) = W (T)^kU(T). 
We calculate the unitary evolution operator U(T) in the 
adiabatic T — > 00 limit in Appendix |Bj by integrating 
over the Berry connection. In the limit U min — > we 
recover the result of Ref . [TS] , 

1k{T) = -s fc Sfc'7fc'(0). (20) 

The result after the three steps is that the Majorana's 
at sites 1 and 2 are switched, with a difference in sign, 

7i(3T) = -sis 2 7 2 (0), 7 2(3T) = s lS27 i(0). (21) 

The corresponding unitary time evolution operator, 

U (3T) = — (l + S1S27172) = exp (^sis 2 7i7 2 ) , (22) 

has the usual form of an adiabatic braiding operation 
[llj . For a nonzero J7 m in the coefficient 7r/4 in the ex- 
ponent acquires corrections of order J7 m i n /J7 max , see Ap- 
pendix [B] 

If one repeats the entire braiding operation, the Majo- 
rana's 1 and 2 have returned to their original positions 
but the final state differs from the initial state by a uni- 
tary operator U(3T) 2 = S1S27172 and not just by a phase 
factor. That is the hallmark of non-Abelian statistics 



entirely different origin (gate-tunable tunnel coupling), 
their Hamiltonian has the same form. We can therefore 
directly adapt their braiding protocol to our control pa- 
rameters. 

We have three fluxes $1, <J> 2 , $3 to control the cou- 
plings. The braiding operation consists of three steps, 
see Table [I] and Fig. [4] (Ref. [TH] had more steps, in- 
volving 6 rather than 4 Majorana's.) At the beginning 
and at the end of each step two of the couplings are off 
(tfrfc = 0) and one coupling is on (\$k\ = ^max)- We de- 
note by Okk' the step of the operation that switches the 
coupling that is on from k to k' . This is done by first 
increasing \§k'\ from to <E> m ax and then decreasing \&k\ 
from <I> max to 0, keeping the third flux fixed at 0. 

During this entire process the degeneracy of the ground 
state remains unchanged (twofold degenerate), which is 
a necessary condition for an adiabatic operation. If, in- 
stead, we would first have first decreased |$fc| and then 
increased |3>fe'|, the ground state degeneracy would have 
switched from two to four at some point during the pro- 
cess, precluding adiabaticity. 

We start from coupling 3 on and couplings 1,2 off. The 



IV. DISCUSSION 

In summary, we have proposed a way to perform non- 
Abelian braiding operations on Majorana fermions, by 
controlling their Coulomb coupling via the magnetic flux 
through a Josephson junction. Majorana fermions are 
themselves charge-neutral particles (because they are 
their own antiparticlc) , so one may ask how there can 
be any Coulomb coupling at all. The answer is that the 
state of a pair of Majorana fermions in a superconducting 
island depends on the parity of the number of electrons 
on that island, and it is this dependence on the electri- 
cal charge modulo 2e which provides an electromagnetic 
handle on the Majorana's. 

The Coulomb coupling can be made exponentially 
small by passing Cooper pairs through a Josephson junc- 
tion between the island and a bulk (grounded) supercon- 
ductor. The control parameter is the flux <E> through the 
junction, so it is purely magnetic. This is a key differ- 
ence with braiding by electrostatically controlled tunnel 
couplings of Majorana fermions [18]. Gate voltages tend 
to be screened quite efficiently by the superconductor, so 
magnetic control is advantageous. Another advantage is 
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that the dependence of the Coulomb coupling on the flux 
is governed by macroscopic electrical properties (capaci- 
tance of the island, resistance of the Josephson junction). 
Tunnel couplings, in contrast, require microscopic input 
(separation of the Majorana fermions on the scale of the 
Fermi wave length), so they tend to be more difficult to 
control. 

Both Ref. [18] and the present proposal share the fea- 
ture that the gap of the topological superconductor is not 
closed during the braiding operation. (The measurement- 
based approach to braiding also falls in this category 
[25] •) Two other proposals [T71QI5] braid the Majorana's 
by inducing a topological phase transition (either by elec- 
trical or by magnetic means) in parts of the system. Since 
the excitation gap closes at the phase transition, this may 
be problematic for the required adiabaticity of the oper- 
ation. 

The braiding operation is called topologically pro- 
tected, because it depends on the off/on sequence of 
the Coulomb couplings, and not on details of the tim- 
ing of the sequence. As in any physical realization of a 
mathematical concept, there are sources of error. Non- 
adiabaticity of the operation is one source of error, stud- 
ied in Ref. |26j . Low- lying sub-gap excitations in the 
superconducting island break the adiabatic evolution by 
transitions which change the fermion parity of the Majo- 
rana's. 

Another source of error, studied in Appendix[Bj is gov- 
erned by the off/on ratio J7 m i n /[/ max of the Coulomb 
coupling. This ratio depends exponentially on the ra- 
tio of the charging energy Eq and the Josephson en- 
ergy Ej of the junction to the bulk superconductor. A 
value Ej/Ec — 50 is not unrealistic 7J, corresponding 
to J7 mi „/[/ ma x - 10~ 5 . 

The sign of the Coulomb coupling in the on state 
can be arbitrary, as long as it does not change during 
the braiding operation. Since J7 max oc cos(7r<7j nc i/e), any 
change in the induced charge by ±e will spoil the opera- 
tion. The time scale for this quasiparticle poisoning can 
be milliseconds [57], so this does not seem to present a 
serious obstacle. 

A universal quantum computation using Majorana 
fermions requires, in addition to braiding, the capabil- 
ities for single-qubit rotation and read-out of up to four 
Majorana's pQ. The combination of Ref. [5] with the 
present proposal provides a scheme for all three opera- 
tions, based on the interface of a topological qubit and 
a superconducting charge qubit. This is not a topolog- 
ical quantum computer, since single-qubit rotations of 
Majorana fermions lack topological protection. But by 
including the topologically protected braiding operations 
one can improve the tolerance for errors of the entire 
computation by orders of magnitude (error rates as large 
as 10% are permitted [5])- 

A sketch of a complete device is shown in Fig. [5j 




FIG. 5: Josephson junction array containing Majorana 
fermions. The magnetic flux through a split Josephson junc- 
tion controls the Coulomb coupling on each superconducting 
island. This device allows one to perform the three types of 
operations on topological qubits needed for a universal quan- 
tum computer: read-out, rotation, and braiding. All opera- 
tions are controlled magnetically, no gate voltages are needed. 
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Appendix A: Derivation of the Majorana- Coulomb 
Hamiltonian 



1. Single island 

Considering first a single island, we start from the 
Cooper pair box Hamiltonian ([2| with the parity con- 
straint ([5]) on the eigenstates. Following Ref. [5T], it is 
convenient to remove the constraint by the unitary trans- 
formation 



H = Q^HCl, n = cxp[z(l -T)(j)/4]. 



(Al) 



The transformed wave function ^((f>) = ^^(cj)) is then 
27r-periodic, without any constraint. The parity operator 
V appears in the transformed Hamiltonian, 



1 



H=—(Q+\e{l-V)+q^ A ) -Ej cos. 



(A2) 



For a single junction the parity is conserved, so eigen- 
states of H are also eigenstates of V and we may treat 
the operator V as a number. Eq. ( A2 1 is therefore the 



Hamiltonian of a Cooper pair box with effective induced 
charge q c s = q- ln d + e (l — 7 ? )/2- The expression for the 
ground state energy in the Josephson regime Ej Eq 



G 



is in the literature [551125] . 

Aground = — Ej + \J 2EqEj 

- W{E c E J /2ir 2 ) 1 / 4 e -^ 8E -' /Ec cos(^ cff /e). (A3) 

The first term —Ej is the minimal Josephson en- 
ergy at (^ m i n = 0. Zero-point motion, with Joseph- 
son plasma frequency u> p — \/8EcEj /h, adds the sec- 
ond term \/2EqCj = ^huj p . The third term is due to 
quantum phase slips with transition amplitudes t± ~ 
exp(±z7r<7 c ff I e) */ hui p Ej exp(— hui p /Ej) by which <fi incre- 
ments by ±27r. 

Using V 2 — f, the ground state energy (A3) may be 
written in the form 



E, 



ground 



= -Ej + x/2E c Ej - UV, 



(A4) 



with U defined in Eq. Q . Higher levels are separated by 
an energy Hui p , which is large compared to U for Ej 
Ec- We may therefore identify -E gr0 und = Heft with the 
effective low-energy Hamiltonian of a single island in the 
large- -Ej limit. 



2. Multiple islands 

We now turn to the case of multiple islands with tun- 
nel coupling. To be definite we take the geometry of 
two islands shown in Fig. [2] The full Hamiltonian is 
H = H\ + H 2 + H12, where Hi and H 2 are two copies of 
the Cooper box Hamiltonian ([2| and Hi 2 is the tunnel 
coupling from Eq. 

To obtain 27r-periodicity in both phases fa and fa, we 
make the unitary transformation H = WHCl with 



Q _ e »(l-Pi)^ I /4 e t(l-7'2)02/4 



(A5) 



The Cooper pair box Hamiltonians are transformed into 

1 2 

H k = 2g (Qk + eq k + q ind ,k) - Ej.k cos fa, (A6) 

with qk = |(1 — ~Pk)- The tunnel coupling transforms 
into 

#12 = ie-^^^iae^ 2 ^ 2 e M*i+»0/a*o + H . c ., (A7) 

where Ti2 = £-£7^712721 and H.c. stands for Hermitian 
conjugate. Since e lq( ^ — cos0 + iqs'mcj), the transformed 
tunnel coupling H12 is 27r-periodic in fa and fa. 



For Ej ^> Ec the phases remain close to the value 
which minimizes the sum of the Josephson energies to the 
bulk superconductor and between the islands. To leading 
order in Em/Ej <C 1 this minimal energy is given by 

£min = - Ej.i - E J<2 + T 12 cos[tt($ 1 + $ 2 )/2$ ] 

+ 0{E\ [ /E J ). (A8) 

The Josephson coupling of the islands changes the plasma 
frequency u) Pl k for phase fa by a factor 1 + 0(Em I Ej), 
so the zero-point motion energy is 

\hw p . k = yj2E c Ej M +E M x OiEc/Ej) 1 ' 2 . (A9) 

The transition amplitudes r± for quantum phase slips of 
phase fa are similarly affected, 

T±, k = -UkVk+EMe-^x/^xOiEc/Ej) 1 ^. (A10) 



These are the contributions to the effective Hamil- 
tonian H cS = £ min + J2k(h^P.k + T +.k + T-,fc) for 
Ej ~> Ec, Em, 



H eS = [-UJ>i - U 2 V 2 + T 12 cos[^($ 1 + $ 2 )/2$ ] 

x [1 + 0{E M /Ej)] + const. (All) 



Eq. ( 12 ) in the main text generalizes this expression for 



two islands to an arbitrary number of coupled islands. 



Appendix B: Calculation of the Berry phase of the 
braiding operation 



We evaluate the unitary evolution operator U of the 
braiding operation in the adiabatic limit. This amounts 
to a calculation of the non-Abelian Berry phase (inte- 
gral of Berry connection) of the cyclic variation of the 
interaction Hamiltonian i?i nt (Ai, A 2 , A3). 

In the Fock basis |00), |01), |10), |11) the interaction 
Hamiltonian ( 18 ) of 4 Majorana fermions is given by 



the occupation number of the two fermionic operators 
c i = (71 — ^72) /2 and c 2 = (70 — 273) /2. It takes the 
form 



Hn 



( -a 3 

A 3 
iAi - A 2 
\ 1A1 - A 2 





7Ai - A 2 
-A 3 




-iAi - A 2 \ 



A3 



(Bl) 
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FIG. 6: The braiding path in three-dimensional parameter 
space along which the Berry phase is evaluated. This path 
corresponds to the flux values in Table[I] with couplings = 
A min for $ fc = and A fc = A max for |$ fc J = $ max . The ratio 
A m i n /A max in the figure is exaggerated for clarity. 



The eigenvalues ar e doubly degenerate at energy ±e = 



A| (up to a flux-dependent offset, which 



only contributes an overall phase factor to the evolu- 
tion operator). The two degenerate ground states at — s 
are distinguished by an even (e) or odd (o) quasiparticle 
number, 



E-Ao 



2s 



£ + A 3 



2c 



/ . £ + A 3 \ 
Ai + iA 2 



V 1 / 

/ \ 

■ £~A 3 

S Ai + iA a 
1 

V / 



(B2a) 



(B2b) 



This parameterization is smooth and continuous except 
along the line Ai = A2 = 0. 

If we avoid this line the Berry connection can be read- 



ily evaluated. It consists of three anti-Hermitian 2x2 
matrices Ak, 











(B3) 



Off-diagonal terms in At are zero because of global parity 
conservation. Explicitly, we have 
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(B5) 



(B6) 



A closed path C in parameter space has Berry phase 



U 



exp ^- j> A k dAi^j 



(B7) 



The path C corresponding to the braiding operation in 
Fig. [1] and Table [I] is shown in Fig. [6] We take all cou- 
plings A/- positive, varying between a minimal value A m ; n 
and maximal value A max . The parametrization (B2| is 
well-defined along the entire contour. 
The contour integral evaluates to 



U = exp 

3 A E 
V2A n 



u 



1 — — e cr 








A, 



(B8) 
(B9) 



The limit A m ; n /A max — > corresponds to the braiding 
operator (22) in the main text (with s±, s 2 > and a z — 
1 - 2c\ci = 17172). 
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